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<~| ' Abstract 

We consider the Euler equation for an incompressible fluid on a three 

dimensional torus, and the construction of its solution as a power series in 

time. We point out some general facts on this subject, from convergence issues 

for the power series to the role of symmetries of the initial datum. We then 

^ . turn the attention to a paper by Behr, Necas and Wu [5]; here, the authors 

chose a very simple Fourier polynomial as an initial datum for the Euler 

equation and analyzed the power series in time for the solution, determining 

VO . the first 35 terms by computer algebra. Their calculations suggested for the 

£f) \ series a finite convergence radius T3 in the i? 3 Sobolev space, with 0.32 < 

T3 < 0.35; they regarded this as an indication that the solution of the Euler 

equation blows up. 

We have repeated the calculations of [5J, using again computer algebra; 

the order has been increased from 35 to 52, using the symmetries of the initial 

rS ■ datum to speed up computations. As for T3, our results agree with the original 

computations of [5] (yielding in fact to conjecture that 0.32 < r 3 < 0.33). 

Moreover, our analysis supports the following conclusions: 

(a) The finiteness of T3 is not at all an indication of a possible blow-up. 

(b) There is a strong indication that the solution of the Euler equation does 
not blow up at a time close to T3. In fact, the solution is likely to exist, at 
least, up to a time #3 > 0.47. 

(c) Pade analysis gives a rather weak indication that the solution might blow 
up at a later time. 

Keywords: Euler equation, existence and regularity theory, blow-up, symbolic 

computation. 
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1 Introduction 

Let us consider the three-dimensional Euler equation for a homogeneous incompress- 
ible fluid (of unit density) with initial datum u , i.e., 

du 

— = —u»vu — Vp , u(x, 0) = uo(x) . 

ot 

The unknown is the divergence free velocity field (x,t) \-¥ u(x,t); we assume peri- 
odic boundary conditions, soi= (x\,X2,Xz) ranges in the three dimensional torus 
(R/27rZ) 3 . In the sequel, we often write u(t) for the function x h- >■ u(x,t). 

One can try a solution of the above Cauchy problem in the form of a power 
series u{t) = Y^ J j=o u j^ (with Uj = Uj(x)); such power series have been the object 
of rather extensive investigations. Morf et al [T5] , Frisch [12], Brachet et al [S], Pelz 
[20J , and other authors (see the bibliography of the cited references) have constructed 
by computer algebra techniques many terms of the power series for specific initial 
data, consisting of simple Fourier polynomials; more precisely, the data analyzed in 
these works are the so-called "Taylor-Green vortex" , and other vortices proposed by 
Kida [H]. The cited authors have also discussed the possibility of a blow-up (i.e., 
finite-time divergence of u(t)) on the grounds of their computer algebra calculations. 
Another initial datum (again a Fourier polynomial) has been considered by Behr, 
Necas and Wu [5] ; these authors have constructed 35 terms of the power series, and 
claimed to have found evidence for a blow-up of the solution; however, in comparison 
with the vortices of Taylor-Green and Kida, the Behr-Necas-Wu initial datum has 
received less attention in the literature. 

The purpose of the present paper is twofold, 
(i) First of all, we wish to point out a number of general facts on the solutions of the 
Euler equation and, in particular, on the convergence of the power series Ylj^o UjP; 
this is the subject of Sections |2] and El Here we report some results extracted from 
the existing literature on the Euler equation in spaces of analytic functions and/or 
in Sobolev spaces; in addition to these results, we present some remarks of ours and 
propose a general treatment to discuss the symmetries of the initial datum and their 
effects on the solution of the Euler equation. We think it is not useless to collect all 
these theoretical statements in a unifying framework, suitable for direct application 
to computer algebra calculations. 

(ii) Our second aim is to reanalyze the power series for the Behr-Necas-Wu initial 
datum, both from the theoretical and from the computational viewpoint; this is the 
subject of Sections HI |5]and[6l First of all we apply to the Behr-Necas-Wu case our 
general setting for the symmetries of the initial datum. We calculate the symmetry 
subgroup of the Behr-Necas-Wu datum (that we recognize to be the dihedral group 
of order 6; this group also determines what we call the pseudo-symmetry space of 
the datum). 



With these premises, we present a novel computation of the power series for the 
Behr-Necas-Wu datum, based on a Python program written for this purpose; this 
computation attains the order 52. The Python program uses an exact representation 
of rational numbers as ratios of integer, so as not to introduce rounding errors; 
furthermore, it employs the symmetries of the initial datum to reduce the amount 
of calculations. 

The results of such computations can be analyzed using the theoretical frame- 
work of Sections |2] and [3j Our conclusions are the following: 

(a) We agree with the estimates of [5], according to which the power series under 
consideration has a convergence radius 0.32 < r 3 < 0.35 in the Sobolev space H 3 ; 
in fact, our computations suggest 0.32 < r 3 < 0.33. However, we disagree from the 
authors of [5] when they interpret the fmiteness of r 3 as indicating a blow-up of the 
solution. 

(b) On the contrary, we give evidence that the solution u(t) of the Euler equation 
exists for t sensibly larger than r 3 . In fact, analyzing the power series for the squared 
Sobolev norm ||w(t)|||, we find a strong indication for a convergence radius 6* 3 such 
that 0.47 < 63 < 0.50. By a general criterion a la Beale-Kato-Majda, this implies 
that the solution of the Euler equation exists, at least, up to time 83. 

The final part of our analysis concerns an alternative approach to estimate 83, 
and the possibility that u(t) blows up at times larger than 83. In connection with this 
problem we use the idea (employed in [8] p2] [15] [20] for different initial data) to 
construct the Pade approximants for the (squared) Sobolev norms and analyze their 
singularities. In particular, we construct the diagonal Pade approximants \p/p](t) 
for ||w(t)|||, up to p = 26. For most of them the complex singularities of minimum 
modulus have modulus ~ 0.5; this fact yields new evidence for the previous estimate 
on 83. Moreover, most of these Pade approximants have real singularities, distributed 
rather erratically; analyzing them in terms of mean value and variance, we obtain a 
somehow weak indication that: 

(c) u(t) might blow up for t — > T - (and t — > (— T) + ), for some T such that 0.56 < 
T < 0.73. 

The blow-up problem can be studied as well in terms of D-log Pade approximants; 
these do not give a clear indication supporting conjecture (c), as briefly explained 
at the end of the paper. In general, much caution is recommended about the Euler 
equation and blow-up predictions via Pade analysis: for example, in the case of the 
Taylor-Green vortex the Pade approximants exhibit real singularities [12] [15], but 
the numerical solution of the Euler equation by spectral methods raises doubts on 
the actual existence of a blow-up [7] [H] . 

Connections with other works. Concluding this Introduction, to put the subject 
of this paper into a wider perspective we wish to mention that there are general 
methods of functional analysis to obtain quantitive lower bounds on the time of 
existence T of the solution of the Euler (or Navier-Stokes) Cauchy problem, from 



the a posteriori analysis of an approximate solution; such lower bounds are certain 
(i.e., non conjectural). 

Derivations of such a posteriori lower bounds have been given in [11] [17] |19| . 
The last of these works gives an algorithm to obtain these lower bounds analyzing 
any approximate solution of the Euler (or Navier-Stokes) Cauchy problem via a 
suitable differential inequality, called therein the " control inequality" . 

Again in [19], a preliminary analysis of the Euler (and Navier-Stokes) equations 
with the Behr-Necas-Wu initial datum has been performed, using for the solution a 
Galerkin approximation with very few Fourier modes. This approximant, combined 
with the control inequality, gives for the Euler equation with this datum a (poor, 
but certain) lower bound T > 0.066 for the time of existence in H 3 (the same 
approach, applied to the Navier-Stokes equations, grants T = +oo when the viscosity 
coefficient is above an explicit threshold). We plan to continue in future works the 
analysis of the Behr-Necas-Wu intial datum, combining the control equation of [19] 
with approximation methods based on extensive automatic computations such as 
the ones presented in this paper. 

2 The Cauchy problem for the Euler equation on 
a torus 

Preliminaries. If a = (a s ),b = (b s ) are elements of R 3 or C 3 , we intend a»b := 
Yls=i a sb s - We indicate with the complex conjugate (and we let it act componen- 



twise on elements of C 3 ); we put \a\ := y/a»a = J^] s=1 \a s \ 2 . 
The Cauchy problem for the incompressible Euler equation is 

— = -tt.Vu-Vp, u(x,0) = Uq{x) , (2.1) 

where: u = u(x,t) is the divergence free velocity field; the space variables x = 

(x s ) 5= i ) 2,3 belong to the torus T 3 := (R/27rZ) 3 ; (u»V«) r := J2 s =i u sd s u r (r = 

1,2,3); p = p(x,t) is the pressure; u = u (x) is the initial datum. As well known, 

the pressure can be eliminated from ( 12.11) using the Leray projection £ onto the 

space of divergence free vector fields; this allows to rewrite the evolution equation 

in (12. ip as du/dt = —£(u»Vu). In this way, we obtain for the Cauchy problem the 

final form 

du 

— = ?(u,u), u(.,0) = u , (2.2) 

where we have written 7 for the bilinear map sending two (sufficiently regular) 
vector fields v,w : T 3 — y R 3 into the vector field 



') 



P(v,w) :=-£(«. Vw) . (2.3) 



In this framework, it is convenient to associate to a vector field v : T 3 — > R 3 the 
Fourier components Vk := (27r) -3 J* T3 dxe~ lh ' x v(x) G C 3 , so that 






v ( x ) = ^ Vk e tk ' x . (2.4) 

fcez 3 

Due to the reality of v , we have V-k = Vk, and v is divergence free iff k»v k = for all 
k. With 7 as above and v,w two vector fields, the Fourier components of y(t>,u;) 
are 

ff(t>, w) k = -i } J Vh»(k - h) ZkWk-h (2.5) 

hez 3 
where 2k '■ C 3 — > C 3 is the projection on the orthogonal complement of k (£*;C := 
c -(k.c)k/\k\ 2 if k ^0; £ c:=c). 

In the above, we have introduced the setting for the Euler equation in an informal 
way; to go on, it is necessary to specify the functional spaces to which the velocity 
fields (at any time) are supposed to belong. 

The expression "a vector field T 3 — > R 3 " can be understood, with very wide gen- 
erality, as "an R 3 - valued distribution on T 3 " (see, e.g., [H]); we write D'(T 3 , R 3 ) = 
D' for the space of such distributions. Any v G -D'(T 3 ,R 3 ) can be differentiated 
in the distributional sense and has a (weakly convergent) Fourier expansion with 
coefficients Vk G C 3 , such that vj: = t>_ fc . 

To construct the full setting for the Euler equation, one must confine the atten- 
tion to much smaller functional spaces of vector fields. For our purposes, two cases 
are important: 

(i) The Sobolev space H n of zero mean, divergence free vector fields of any order 
n G [0, +oo). This is defined in terms of the space L 2 (T 3 ,R 3 ) = L 2 of square 
integrable vector fields v : T 3 — > R 3 , equipped with the inner product (v \w)l,2 := 



(2%) 3 J T 3 v(x)»w(x)dx and with the induced norm ||v||l 2 = (27r) 3 ^ 2 -\// T 3 \v(x)\ 2 dx 

(note the term (27r) -3 in the inner product, used systematically in the sequel). By 
definition, 

F n (T 3 , R 3 ) = H n := iv E D' | V^A n v e L 2 , f v dx = 0, divv = o\ (2.6) 

= iv e D' | ^2 \ k \ 2n \ v k\ 2 < +oo, ^o = 0, k.v k = 0J . 
fcez 3 

(In the above a/— A indicates the power of order n/2 of minus the Laplacian; by 
definition (y/-A n v) k = \k\ n v k for each v G D'. Note that H n C L 2 for all n > 0.) 
H n is a Hilbert space with the inner product 

(v\w) n := (V^A n v\V^X l w) L 2 = Y, \k\ 2n Vk-'W k , (2.7) 

fcez 3 



inducing the norm 



IMU = \\%f=K n v\\ I? = Y. W 2n M 2 • ( 2 - 8 ) 

It is known that 7 sends continuously H n x H n+1 into H n , for all n G (3/2, +00). 
(ii) The space of C^ (i.e., analytic) zero mean, divergence free vector fields on T 3 ; 
this is 

A(T 3 , R 3 ) =A:= Ive C W (T 3 , R 3 ) | f v dx = 0, divw = } (2.9) 

v € D' I liminf |t; fe |~T*TT+l*£l+M > 1, y = 0, k.v k = 
fcez 3 , fc-s-oo 

1 

(intending Fn+Rl+ra := -)-oo. The Fourier representation in f |2.9j) mimics the 

description of analytic functions on the torus in [16] , which is also a useful reference 

for what follows). One has 

A = U pG (i i+00 ).Ap , (2-10) 

A p := \v G D' I liminf \v k \~ i fc ii+i fc 2i+i fc 3i > p, Vo = 0, k,v k = o) ; 
I fcez 3 ,fe-i>oo J 

each yi p is a vector subspace of A Let us introduce the annulus K p := {z G 
C I 1/p ^ |z| ^ p} and its power K p := {z = (zi,Z2,z 3 ) G C 3 | Zi,Z2,Z3 G K p }. 
For t> G yip, the series ^ fcGZ 3 UfcZ fe converges in C 3 for each z G K z p (intending 
z k := z 1 1 z 2 2 z 3 3 ); the function 2; 1— > ^2k&z 3VkZ,k ls holomorphic on the inner part of 
K p and continuous on K p , so we can define 



\ P ■= sup 

zGA-3 



J2 Vkzk 



kez 3 



(2.11) 



HI || I p is a norm on A p and makes it a Banach space. One equips A with the inductive 
limit topology of the collection of Banach spaces {(A p , ||| ||| p ) | p G (1, +00)}: this is 
the finest locally convex topology on A making continuous each embedding A p ^-> A. 
(Besides [16] . see [23] for the general theory of inductive limits.) A is continuously 
embedded into each Sobolev space H n ; the map CP is continuous from A x A to A. 

Basic results on local existence and uniqueness. We start from the Sobolev 
framework, choosing 

ne (5/2, +00). (2.12) 

In the sequel, an if "-solution of the Euler equation, or of the Euler Cauchy problem, 
means a map 

ueC((-7,T),H n )nC 1 ((~7,T),H n - 1 ) (2.13) 



(T, T G (0, +00]) fulfilling the Euler equation, or its Cauchy problem with a suitable 
initial condition u(0) = Uq. The following statement is well known: 

2.1 Proposition. For n as above and any initial datum uq G H n , the following 
holds. 

(i) The Cauchy problem Ii2. || ) has a unique maximal (i.e., not extendable) H n - 
solution u of domain (—7, T), for suitable T = T(u ), 7 = 7(u ) G (0, +00]. 
(ii) (Beale-Kato-Majda criterion, Sobolev version). If T < +00, one has 

T 

dt\\u(t)\\ n = +00 , (2.14) 



a fact implying 

limsup ||w(t)|| n = +00 . (2-15) 

Similar results hold if7< +00 ; considering the integral from —7 to and the limit 
fort-* (-T)+. 

Proof, (i) See @] [13]. 

(ii) See [3]. Indeed, here it is shown that T < +00 implies J Q dt\\rotu(t)\\L°° = 
+00; however, ||rotw(t)||i«, ^ const. ||«(t)|| n by the Sobolev imbedding inequalities, 
whence Eq. fj2.14j) and its obvious consequence (12.151) . The behavior of u at time 
— 7 is analyzed similarly. □ 

If T < +00, the solution u is said to blow up at time T. Similarly, if 7 < +00 
we say that u blows up at — T. Many statements presented in the sequel on the 
possibility of blow-up at T have obvious reformulations regarding — 7. 

2.2 Remark. The Beale-Kato-Majda criterion (12.141) yields the following state- 
ment, in case of blow-up with a power law: 

if |Kt)|| n ~ 7^ ^ for t -> T~ (with U,a> 0), then a ^ 1 . (2.16) 

[1 — t) 

In the case of the Euler equation on R 3 , it was recently shown in [10] . Theorem 1.3 
that the blow-up at T implies the following, for any n > 5/2: 

IN*)lln> (T _^ )2 n/5 fo^ close to T, U = U n (\\u \\ L2 ). (2.17) 

This estimate might hold as well for the framework of the present paper, i.e., for 
the Euler equation on the torus T 3 (however, the extendability of ( I2.17P to T 3 is 
immaterial for the purposes of this paper). □ 



Let us pass to the C u (= analytic) framework; what follows assumes some general 
notions from the theory of analytic functions from R to locally convex spaces, for 
which we refer to [6] §3. Let A be the space f 12 . Q j) ; in the sequel, an .A-solution of 
the Euler equation, or of the Euler Cauchy problem, means a map 

ueC"{{-7,T),A) (2.18) 

(7, T G (0, +oo]) fulfilling the Euler equation, or its Cauchy problem with a suitable 
initial condition u(0) = uq. Let us report a known result. 

2.3 Proposition. For any initial datum uq G A, the following holds, 
(i) Problem \2. || ) has a unique maximal (i.e., non extendable) A-solution of domain 
(-T,T), for suitable T = T(u ),7 = 7(u ) G (0,+oo]. 

(ii) For any n G (5/2, +oo), this coincides with the maximal H n -solution of the 
Cauchy problem with the same datum (and thus, if T < +oo ; it fulfills Eqs. \2.1J$ 
l\2.15\) ; a similar result holds if7< +00^. 

Proof, (i) See [2J, Theorem III. 2, page 264 (this is a result of existence and uni- 
queness on sufficiently small time intervals, from which one infers via standard 
arguments existence and uniqueness of the maximal solution), 
(ii) See [3], especially Remark 2.1, page 414. □ 

Assuming again uq G A, and choosing any n G [0, +00), we conclude with two 

remarks. 

(i) By the continuous embedding of A into H n , the function u of the last proposition 

isalsoinC w ((-T,T),F n ). 

(ii) Consider the function 

(-T,T)->R, t^\\u(t)\\ 2 n . (2.19) 

This is in C UJ ((—7, T), R), being the composition of the analytic function u : (— 7, T) 
— >■ H n with the continuous quadratic function || ||^ : H n — > R. 

Symmetries of the Euler equation. Let us consider the octahedral group Oh, 
formed by the orthogonal 3x3 matrices with integer entries: 

O h ■= {S G Mat(3 x 3, Z) I S T S = 1 3 } . (2.20) 

In fact, the entries of any such matrix have —1,0 and 1 as the only possible values; 
furthermore, a 3 x 3 matrix S belongs to Oh if and only if 

5 = diag( ei ,e 2 ,e3)Q(a) (2.21) 

e s G {±1} (s = 1, 2, 3) ; Q(o~) the matrix of the permutation a : {1, 2, 3} — > {1, 2, 3}; 

more precisely, Q(a) is the matrix such that (Q(a)c) s = c CT ( s ) for all c G C 3 , s G 



{1, 2, 3}. There are 2 3 = 8 possible choices for the signs e, and 3! = 6 choices for a, 
so Oh has 8 x 6 = 48 elements. Clearly, each S G Oh sends Z 3 into itself. 

To go on, let us denote with Oh x T 3 the Cartesian product Oj, x T 3 , viewed as 
a group with the composition law defined by (p) 

(S, a)(U, b) := (SU, a + Sb) (S,U £ O h ; a,b G T 3 ) . (2.22) 

Of course, the unit of this group is (1,0) (with 1 the identity 3x3 matrix); the 
inverse of a pair (S, a) is (S, a)" 1 = (S T , —S T a). To any element (S, a) of Oh x T 3 is 
associated a "rototranslation" 

E(S,a) :T 3 ^T 3 , x i-4 £(S, a)(x) := Sx + a , (2.23) 

and one checks that the mapping (S,a) H- E(S,a) is a group homomorphism be- 
tween Oh x T 3 and the group of diffeomorphisms of T 3 into itself (with the usual 
composition). 

Now, we take a vector field v in H n (or in A) and an element (S, a) of the group 
Oh x T 3 . We can construct the push-forward £*(S, a)v of v along the mapping 
£(S, a); this is the vector field in H n (or in A), given by 

£*(S,a)v:T 3 ^R 3 , x \-> (£*(S, a) v)(x) = Sv(S T (x - a)) . (2.24) 

One easily checks that Eq. (12.24p actually defines a vector field in H n (or in A), 
with Fourier components 

(£*(S, a) v) k = e~ ia ' k Sv S T k (k G Z 3 ) . (2.25) 

Let us write £*(S f , a) for the map v G H n \-t £*(S, a)v; this is a linear map of H n into 
itself, preserving the inner product ( | ) n , so it is in the group 0(H n ) of orthogonal 
operators of the Hilbert space H n into itself. The mapping 

£,:O ft KT 3 4 0{H n ) , (3, a) i-> £,(5, a) (2.26) 

is a injective group homomorphism, i.e., a faithful orthogonal representation of the 
group Oh x T 3 on the real Hilbert space H n . Alternatively, let us write £*(S", a) for 
the map v G A i— > £*(S", a)t>; this is in the space Iso(.A) of linear and topological 
isomorphisms of A into itself. The map 

£,:OikT 3 4 Iso(yi) , (S, a) i-> £,(5, a) (2.27) 

is an injective group homomorphism, i.e., a faithful linear representation of the group 
Oh x T 3 on the topological vector space A. 



2 This is the semidirect product of the groups Oh and T 3 with respect to the natural homomor- 
phism Oh — > Aut(T 3 ) sending S € Oh into the map b i-> Sb, an automorphism of T 3 . 



Let us relate the previous constructions to the bilinear map CP of the Euler 
equation. From the Fourier representations (12.51) (j2.25p . one easily infers 

?(E*(S, a) v, £*(S, a) w) = E*(S, a) V(y, w) (2.28) 

for all v G H n , w G H n+1 with n > 3/2 (and, in particular, for all v , w G A). Let us 

outline the implications of H2.28[) about the solutions of the Euler equation. In the 

rest of the paragraph, the term "solution" either means an if "-solution (n > 5/2) or 

an .A-solution, and the initial datum uq is chosen consistently in H n or in A. From 

(I2.28P one infers the following, for each (S, a) G Oh x T 3 : 

(i) If u : t G (— T, T) I—?- u(t) is a solution of the Euler equation, we have two more 

solutions 

£*(£, a)u: tE (-7, T) H- £,(5, a)u(t), (2.29) 

- £,(5, a)u(-) : t G (-T, 7) ^ -£*(S, a)u(-t). (2.30) 

(ii) If u : t G (— T, T) i-> w(t) is the maximal solution of the Euler Cauchy problem 
with datum m , then £*(S, o)u is the maximal solution with datum £*(5, o)«o and 
— £*(5, a)w(— •) is the maximal solution with datum — £*(S, a)«o- 
(iii) Let us denote again with u : t G (— T, T) i— >■ w(t) the maximal solution of the 
Cauchy problem with datum u . Then, 

e*(S,a)u = u => £.*(S,a)u(t) = u{t) for tG(-T,T). (2.31) 

- £*(£, a)u = u => 7 = T, - £*(S,a)u{-t) = u{t) for t E (-T,T). (2.32) 

The verification of statements (i) (ii) is straightforward. After this, the implication 
(j2.3ip in (iii) follows noting that £*(5 r , a)u and u are maximal solutions of the 
Cauchy problem with the same datum £*(£, o)uq = u . Similarly, the implication 
(I2.32p follows noting that —E*(S, a)u( — ) and u are maximal solutions of the Cauchy 
problem with the same datum — £*(5', a)u = u . 

Considering the maximal solution u for a datum w i n H n (n > 5/2), and recalling 
that any transformation £*(S f , a) preserves the H n norm, we also obtain from (j2.32p 
the following: 

~E^(S,a)u = u =}► T = T, \\u(-t)\\ n = \\u(t)\\ n for tG(-T,T). (2.33) 

The results in (iii) suggest to consider, for a given datum uq in H n or A, the 
symmetry subgroup 

H(u Q ) := {{S, a) G O h k T 3 | £,(5, a)u = u } (2.34) 

and the pseudo-symmetry space 

%-{u Q ) := {(S, a) G O h x T 3 | - £*(S, a)u = u } (2.35) 

9 



(the first one, being a subgroup of Oj, k T 3 , contains at least the identity element 

(1,0); the second one might be the empty set. The term "isotropy group", often 

employed in place of "symmetry group", will not be used in this paper). 

Let us consider the maximal solution u of the Cauchy problem with a datum w 

(contained in H n for some n > 5/2); from Eqs. fl2.3ip - fl2.33p . we readily obtain the 

following: 

£,(5, a)u{t) = u(t) for all (S, a) G U{u Q ),t G (-T, T) ; (2.36) 

7T(uo)^0 => 7 = T, -Z*(S,a)u(-t)=u(t), ||u(-*)|| n = \\u(t)\\ n (2.37) 

for (S, a) eU-(u ),te (-T,T) . 

For future use, let us introduce the reduced symmetry subgroup and the reduced 
pseudo-symmetry space of the datum u , which are 

7~t-R(uo) := {S G Oh | £*(£, a)uo = Uq for some a G T 3 } , (2.38) 

Hr{uq) := {S G O h | - £*(5, a)u = u for some a G T 3 } . (2.39) 

Let us observe that the set theoretical unions %{uq) WH~(uq) and T-Lr(uq) WH^Uq) 
are subgroups of Oh x T 3 and Oh, respectively. 

As a final remark, useful for the sequel, let us consider the pair (—1,0) G Oh x T 3 , 
noting that £(—1,0) is the space reflection: £(—1,0) (re) = — x for all x G T 3 . One 
easily checks that 

(-l,0)er(uo) ^ K-(uo) = H(uo)(-l,0) = {(-5 , ,o)|(5 , ,o)eH(u )} (2.40) 

(where 7i(u )(— 1,0) stands for the set {(S,a)(— 1,0) | (S, a) G "H(u )}; the last 
equality rests on the identity (S, a)(— 1, 0) = (— S, a)). 

3 Power series in time for the Euler Cauchy prob- 
lem 

Throughout this section, we consider the Euler Cauchy problem with initial datum 
Uq G A. 

Setting up a power series for the solution. Let us try to build the solution of 
the Euler Cauchy problem as a power series 

oo 
3=0 

with coefficients Uj G A, whose convergence has to be discussed later. The zero order 
term in this expansion is the initial datum uq; to determine the other coefficients 

10 



Uj G A, it suffices to substitute the expansion (13.11) into the Euler equation 

and to require equality of the coefficients of the same powers of t in both sides: in 

this way, one easily obtains the recurrence relation 



u j 



3-1 
J 



1 ' 

T ^%^-i) (j = 1,2,3,...) . (3.2) 



e=o 



When applying this recurrence relation for the u/s it can be useful to represent the 
bilinear map CP in terms of Fourier coefficients, as in Eq. ( 12. 5ft . This is especially 
useful if the initial datum uq is a Fourier polynomial, i.e., if uok 7^ only for finitely 
many modes k. In this case, all the iterates Uj (j = 1,2,3,..) are as well Fourier 
polynomials, and the implementation of ( 13. 2 p via the Fourier representation (12. 5p 
always involves sums over finitely many modes. 

In the next section, a large part of our attention will be devoted (for a specific 
datum M ) to the partial sums 



■u 



N 

' v )(;):=$>^, (3.3) 

3=0 



(N = 0, 1, 2, ...) and to the (squared) Sobolev norms 

' (A °(«= E \k\ 2n Wi N \t)\ 2 ■ (3-4) 



\ u 

fcez 3 



Symmetry considerations. Let us consider the symmetry subgroup "H(uo) or 
the pseudo-symmetry space 'H~(uq), see Eqs. (12.341) ( J2.35D . Using the recursive 
definition ( 13.21) of Uj with the invariance property ( 12. 28ft of 7, one easily checks the 
following, for any j G {0, 1, 2, ...}: 

£*(S,a)Uj = Uj for all (S,a) G H(u Q ); (3.5) 

— £*(S,a)uj = (—iyuj for all (S,a) G 'H~{uq). (3.6) 

Of course, the last two equations imply the following, for all N G {0, 1,2,...}, t G R 
and n G [0, +oo): 

£*(S, a)iiN(t) = u N (t) for (S,a) G H(u ); (3.7) 

— £*(S,a)uN(t) = u N (—t) for (S, a) G H _ (u ); (3.8) 

\\u N {t)\\n= \\u N (-t)\\n if«-(«o)^0 (3.9) 

(Eq. (13. 9 p is a consequence of Eq. (13.81) and of the invariance of || || n under the 
transformation &*(S, a)). 
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Due to the Fourier representation (I2.25P for E*(S,a), the equality ( 13 .5p reads 
e -ia»k Suj,s*k — u j,k or, equivalently, 

Uj, S k = e~ ia ' Sk Su j:k for k G Z 3 , (S, a) G H(u ) ; (3.10) 

similarly, Eq. (13 .6p is equivalent to the statement 

u hS k = {-l) 3+1 e- ta ' Sk Su hk for k G Z 3 , (S, a) G ?T(uo) • (3.11) 

In typical applications of the recursion scheme (13. 2p . where uq is a Fourier polynomial 
as well as its iterates Uj, Eqs. (I3.10p (13. lip can be used to speed up the computation 
of the Fourier components of the u/s; in fact, at any given order j, after comput- 
ing a Fourier component Uj^ we immediately obtain from the cited equations the 
components Uj^k for all S in the reduced subgroup or subspace Hr(u ), "H^(tto)- 

Convergence of the power series in A,. From now on, we intend 

r := convergence radius of the series Y^T=o u o^ i n -^ • (3-12) 

Furthermore, 

u : t G (— T, T) h-> u(t) is the maximal .A-solution of the Cauchy problem (3.13) 
(recall that, for any n > 5/2, u is also the maximal H "-solution) . We note that 

oo 

< t ^ 7 A T , u(t) = Y^ u j tJ in A for t G (-r, r) (3.14) 

(with A indicating the minimum). In fact: being analytic, u admits a power series 
representation in a neighborhood of zero; this necessarily coincides with the series 
(13. ip . whose convergence radius r is thus nonzero and fulfills (—t,t) C (— T, T). 

Convergence of the power series in H n . After fixing n G [0, +oo), let us discuss 
the series (13. ip in the Sobolev space H n . To this purpose, we put 



r n := convergence radius of the series Y^jLo u jt 3 m H n ; (3.15) 

the root test gives 

T n = liminf ||«./||„ •* (3.16) 

(intending -1 /- 7 := +oo). With r, T, T, w as before, we claim that 

oo 

t ^ r„ and u{t) = ^ u j fj in H ", for t G (~ T A T n, T A r n ) (3.17) 

3=0 
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(where —7Ar n is the opposite of the minimum 7 A r n ). In fact: the series YlJLo u j^ 
converges to u(t) in A, for t G (— T, r); by the continuous embedding A M- H n this 
series converges to u(t) in if n as well, at least for t G (— r, r); thus r ^ r n . Moreover 
the functions u : (— T, T) — > H n and t G (— r n , r n ) i— )■ YlT=o u jV G if n are analytic 
and coincide on (— r, r); so, by the analytic continuation principle, these functions 
coincide on the intersection of their domains which is (— T A r n , T A r n ). Let as add 
a stronger claim: 

5 °° 

if n > -, t < r„ ^ T A T and w(£) = ^ w^' in # n , for t G (-r„, r„) . (3.18) 

3=0 

In fact, the function t G (-r n , r n ) M- ]T~ Q u i t3 is in C((- r n 5 T n), H n ) n C 1 ^--^, r n ), 

if n_1 ) and solves the Euler Cauchy problem, so it is a restriction of the maximal 

if "-solution, which is u of domain (— T, T); this gives the relations r n ^ T A T and 

M W = Y.f=o u j tJ in # n > for * e (-Tn,T n ). 

Power series for the Sobolev norms of the solution. Let us choose n G 

[0, +oo). The squared norm || ^ji=o M i^ lln = (Sj=o u j'^' I Sj=o u jV)n has the formal 

expansion 



j 



|| ^J M i i J ||£ = ^ z/ ni t J , z/ nj := y^(^|tt 3 -_^) n G R ; (3.19) 

i=o j=0 £=0 

for future use we remark that (p) 

7T(uo) ^ => i/ fti = for all j odd. (3.20) 

Independently of any assumption on H~(uo), let us define 

9 n := convergence radius of the series Yl^Lo Vn ft = lim inf \v n j\ ■ (3.21) 

Let us relate these objects to the convergence radius r n in (I3.15p . to the solution 
u G C^{{— 7, T),A) and to its squared H n norm. We claim that 

+oo 

r n ^ 6 n and \\u(t) \\l = Y^ "njt j for t G (-T A 6 n , T A n ) (3.22) 

3=0 



3 Let us propose a proof of (13.20[) . based directly on the definition (|3.19[) of v n j- If 
W~ (uq) has at least one element (S, a), from Q3.6J) and from the invariance of ( | ) n un- 
der any transformation £.*(S 7 a) we obtain that, for each £ £ {0,...,j}, (ui\iij-e) n = 
((—!)* £*(S, a)u£\(-iy~ e E^(S, a)uj-t} n = (—!)■* (t^|u,-_*)„, whence v nj = {-l) j v ny If j is odd, 



this means u n j = 0. 
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(with — T A 8 n the opposite of 7 A 9 n ). In fact: the expansion u(t) = ^2^^,UjP, 
converging in H n for t G (— r„,r n ), implies r n ^ 6» n and ||u(£)||£ = J2^ u njt j for 
£ G (— T n ,T n ). Moreover the functions t G (— 7,T) H- ||u(£)||£ and t G {-6 n ,6 n ) t-> 
J2j^o u njV are analytic and coincide on (— r n ,r n ), so they coincide everywhere on 
the intersections of their domains, which is (— 7 A 8 n , T A 9 n ). We now add to ( I3.22J) 
a stronger claim: 

e +°° 

ifn>-, r^^^TATand ||u(£)|| r 2 t = ^ z/ nj £ J for £ G (-# n A). (3.23) 

3=0 

Let us prove this claim, assuming for example that 7 AT = T. If it were T < 9 n 
we would infer lim t _> T - ||u(£)||£ = hm^ T - £]+~ z/ nj t J = X)t^v n /P < +oo (the 
first equality would hold due to H3.22|) and T /\6 n = T; the subsequent two relations 
would hold because T would be inside the convergence interval of the series). On the 
other hand, since n > 5/2, the conclusion that lim^^- ||w(£)|| 2 exists finite would 
contradict (j2.15p . 

4 Power series for the Euler equation in a paper 
of Behr, Necas and Wu 

In the paper jS] mentioned above, the authors considered the power series (13. ip for 
the Euler equation on T 3 , with an initial datum u G A given by 

u (x) = Yl u ^ k ' X ' ( 41 ) 

fc=±a,±6,±c 

a := (1,1,0), 6:= (1,0,1), c := (0, 1, 1) ; 
u ,± a := (1, -1, 0) , u ,±b := (1, 0, -1) , u 0)±c := (0, 1, -1) . 

Like «o, all the subsequent terms Uj are Fourier polynomials with rational coefficients 
(□)• Using rules equivalent to (j3.2p f 12 . 5 1) . the terms Wj were determined in [5] by 
computer algebra, for j = 1,2, ...,35. Computations were done with Mathematica 
for j = 1,...,10, and with a C++ program for j = 11,. ..,35 (in the later case, 
approximating the rational coefficients with finite precision decimal numbers). After 
determining the Uj's, the authors fixed their attention on the partial sums 

N 

U W(*):=J>/, 

j=0 



4 For a more precise statement on these cofficients see our discussion of the datum w in the 
next section and, in particular, Eq. ([52 
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whose iV —¥ +00 limit gives the solution u(t) of the Euler Cauchy problem, for all 
t such that the series converges. The previously mentioned computation of the u/s 
made available these partial sums for N = 0,1, ...., 35; the authors of [5] computed 
the (squared) Sobolev norm 

h w ®\\l = E l*l 6 H?°(t)l a 

fcez 3 

for the above values of N, and several values of t. Their main results were the 

following: 

(i) Setting t = 0.32, and analyzing the behavior of ||ti^ iV ' 1 (0.32) || 3 for N from to 

35, the authors found evidence that ||ti^ 7V - ) (0.32)||3 should approach a finite limit for 

JV->- +00. 

(ii) Setting t = 0.35, the authors observed a rapid growth of ||m^(0.35)||3 for iV 

ranging from to 35, a fact suggesting that lini7v_j. +00 ||m^(0.35)||3 = +00. 

(iii) A behavior as in (ii) was found to occur for slightly higher values of t (even 

though the authors suspected some rounding error to appear for t > 0.35). 

The above results suggest that the series ^2^UjP has a finite convergence 
radius r 3 in H 3 , with r 3 e (0.32, 0.35). 

Let us discuss this outcome from the viewpoint of the present paper, denoting 
with u the maximal ^l-solution of the Cauchy problem with this datum and recalling 
that this coincides with the maximal H 3 -solution. The datum uq possesses pseudo- 
symmetries (to be described in the next section); therefore, u has a time symmetric 
domain (—T,T) (in [5] this fact was not explicitly declared, but probably regarded 
as self-evident). According to our Eq. (|3.18p . it is 

r 3 ^T; (4.2) 

in principle, it could be T = +00. In spite of this, the authors of [5] spoke of a 
blow-up at r 3 . 

In the next two sections we present our computations on the power series for the 
Behr-Necas-Wu initial datum, with our interpretation of the results. Even though 
these calculations confirm the "experimental" outcomes (i)-(iii) of [5], we give ev- 
idence that the solution u of the Euler equation does not blow up close to r 3 ; on 
the contrary, computing the power series for ||w(£)|| 3 up the available order we ob- 
tain strong evidence that such a power series has a convergence radius 9 3 such 
that 0.47 < 63 < 0.50, which implies for the time T of existence of u the bound 
T ^ #3 > 0.47. By a subsequent analysis relying on the technique of the Pade ap- 
proximants, we show that a blow-up of u(t) might happen at a time larger than 
0.48: more precisely, these computations give a somehow weak indication that T 
might be finite, with 0.56 < T < 0.73. 
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5 Our approach to the power series of Behr, Necas 
and Wu 

Let us denote again with uq the datum (14.11) and consider its iterates Uj (j = 
1,2, ...), with the corresponding power series; like Uq, all the iterates Uj are Fourier 
polynomials with rational coefficients. Throughout the section, u is the maximal 
.A-solution of the Euler equation with datum uq. 

A closer analysis of the Behr-Necas-Wu initial datum: symmetry prop- 
erties. The symmetry group 7i(u ) and the pseudo-symmetry space H~(u ) (Eqs. 
(I2.34p (I2.35P ) can be explicitly computed. For the first one, we find 

H(u ) = {(1, 0), (1, z 2 ), (A, o x ), (A, a 2 ), (B, a x ), (B, a 2 ), (5.1) 

(C, ci), (C, c 2 ), (D, d), p, c 2 ), (E, 0), (E, i 2 )} 

where 1 is the 3x3 identity matrix, and 



(5.2) 



C := 

furthermore, Z2,Oi, etc., are the following elements of T 3 : 

z 2 := (7r,7r,7r) , oi := (0,0, 7r) , a 2 := (vr, vr, 0) , (5.3) 

ci := (-7T, 0, 0) , c 2 := (0,7T,7r) 

(of course, in the above it is short for it mod. 2ttZ). Let us fix the attention on the 
reduced symmetry subgroup 1-Lr{uo) = {1, A, B, C, D, E}; it is readily checked that 

A 3 = 1,B 2 = 1,(BA) 2 = 1 (5.4) 

C = A 2 , D = AB , E = A 2 B . 

So, 'Hr{uq) has two generators A, B\ the first line in (j5.4p gives a presentation of this 
group in terms of generators and relations, while the second line expresses the other 
elements in terms of A, B. Using Eq. (|5.4j) . one recognizes a group isomorphism 

Hr(uo) ~ D 3 (5.5) 



16 




where the right-hand side indicates the dihedral group of order 3, formed by the 

symmetries of an equilateral triangle p). 

Now we consider the full group H(uo) (with the product ( 12.221) ). It is easy to check 

that 

(A ai f = (1, 0) , (B, ai f = (1, 0) , ((£?, ai )(A, a,)) 2 = (1, 0) , (5.6) 

(A a,) 2 = (C, c 2 ) , (A, a,) 3 = (1, i 2 ) , (A, ai ) 4 = (A, a 2 ) , (A, a,f = {C, Cl ) , 

(A, oi)(S, oi) = (D, c 2 ) , (A ai) 2 (5, oi) = (E, z 2 ) , (A, ax) 3 (B, a x ) = (B, a 2 ) , 

(A ax) 4 ^, oi) = {D, ci) , (A ai) 5 (S, ai) = (£, 0) . 

So, H(uo) has two generators (A °i) an d (-B,ai); the first line in (I5.6P gives a 
presentation of this group in terms of generators and relations, and the subsequent 
lines express the other elements in terms of the generators. One recognizes a group 
isomorphism 

H(uo) ~ D 6 (5.7) 

where the right-hand side indicates the dihedral group of order 6, formed by the 
symmetries of a hexagon (see the previous footnote). 

Let us pass to the pseudo-symmetry space T-L~{uq). One readily checks that this 
contains (— 1, 0) (inducing the space reflection £(—1,0) : x G T 3 i— > —x). From here 
and from the general result ( 12.4=0(1 . one obtains 

K-(u ) = H(u )(-1, 0) = {(-1, 0), (-1, z 2 ), (-A a x ), {-A, o 2 ), (5.8) 

(-B, ai ), (-B, a 2 ), (-C, Ci), (-C, c 2 ), (-£>, Cl ), (-D, c 2 ), (-B, 0), (-E, i 2 )} , 

with A, B, ... and z 2 , a±, ... as in Eqs. (15. 2p (15.31) . 

Some consequences of the previous symmetry results, (i) What we have 
stated in Section [3] for an arbitrary initial datum holds, in particular, for the present 
datum m : the symmetries or pseudo-symmetries of uq can be used to speed up the 
computation of the Fourier components of any iterate Uj. More precisely, if we know 
the Fourier component Uj t k for some k, using Eqs. ( 13.101) ( 13.111) we readily obtain 
the components Uj^k f° r ai l S G 'Hr{uq) U 7i~^(u Q ). 

(ii) As already noted, the pseudo-symmetry space - H _ (m ) contains (—1,0), corre- 
sponding to the space reflection. In terms of Fourier coefficients, the relation (13. lip 
with (S, a) = (—1, 0) takes the form Uj-k — (~^-Y u j,k f° r j — 0, 1, 2, ... and k G Z 3 . 
On the other hand, any iterate Uj is a real vector field, thus Uj_k = uj^k] in con- 



5 For any integer n £ {3,4, ....}, one denotes with D n the dihedral group of order n; this is 
formed by the orthogonal transformations of the Euclidean plane R 2 into itself which preserve a 
regular polygon with n sides, centered at the origin. Denoting with id the identity map, with a 
the rotation of an angle 2tt/ti and with b the reflection about anyone of the n symmetry axes of 
the polygon, one finds that a, b are generators of D n and fulfill the relations a n — id, b 2 = id, 
(ba) 2 = id. The elements of D„ are 2n, and coincide with id, a, a 2 , ...,a™A b, ab, a 2 b, ..., a n ~ 1 b. 
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elusion ujj: = (—lyujk, which indicates that Uj t k is real for j even, and imaginary 
for j odd. Taking into account that the coefficients Uj t k are rational in any case, we 
conclude the following for each fcgZ 3 : 

u jt k e Q 3 for j = 0,2,4, ... ; u jtk G iQ 3 for j = 1, 3, 5, ... . (5.9) 

(iii) In the sequel we are often interested in the partial sums u^ N \t) := J2j=o u j^ 
and in their norms ||tt^(t)|| n , especially for n — 3. Since H~(u ) ^ 0, as in (J3.9P 
we have ||wjv(*)||n = ||«Ar(-*)|U- 

(iv) Independently of any convergence consideration about the power series ^ 7 -^ UjP . 
the result 7i~(u ) ^ also ensures that the (maximal A-) solution u of the Euler 
equation with datum uq has a symmetric domain (—T,T) (recall Eq. (I2.32p ). 

Describing our computations. We have considered again the power series (13.1 j) 
for the datum uo', to deal with this series we have written a program in Python, 
using the package gmpy [23] for fast arithmetics on rational numbers. This program 
implements Eq. (J2.5P for CP and the recursion rule (J3.2P ; moreover, it takes into 
account the dihedral symmetries (and pseudosymmetries) of u to speed up compu- 
tations. The program has been run to compute the terms Uj for j = 1, ...,52 (Jj). 
Calculations have been performed on a PC with an Intel Core i7 CPU 860 at 2.8GHz 
and an 8GB RAM. The CPU time for Uj has been, for example: 1 second for j = 10, 
one minute for j = 20, half an hour for j = 30, 7 hours for j = 40 and 85 hours for 
j = 52. Differently from [5], for all orders up to j = 52 the Fourier coefficients Ujk 
of Uj have been represented as elements of Q 3 or zQ 3 ; so, no rounding errors related 
to finite precision arithmetics have been introduced in the calculation of the power 
series. 

From the Uj's one determines the squared norms ||mj||| = Xlfcez 3 l^| 6 l M i,fc| 2 ; ^ ne 
partial sums u^ N \t) := J2j=o u j^ an d their squared norms ||w^(i)ll3 = Sfcez 3 
l^| 6 l w fe Wl 2 (^ = !)•••) 52). Each ||Wj||| is a rational number and H'M^^Ha is a 
polynomial of order 2N in t, with rational coefficients, containing only even powers 
of t; furthermore, the coefficients of t° and t 2N in ||w^(t)||§ are ||wo|li an d ll M v|l3) 
respectively. 

Our computations of the above norms, up to j = 52 or N = 52, have been done 
using the previously mentioned Python program. These calculations have been 
relatively quick: for example, the computation of ||w^(£)||3 has required a CPU 



6 To test the reliability of this program, the calculation of some of the u's has been checked 
in two independent ways. These checks have been done by means of other two programs, which 
implement Eqs. (|2.5|) p.2[) accepting as an initial datum uq any Fourier polynomial; these do not 
refer to any symmetry property of uq. The first of these programs, written in Mathematica, has 
been used to compute the Uj's up to order j — 13; the second program, written in Python, has 
been used for a calculation up to j = 43. 
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time of about 3 hours. As first examples of our results, we report the following ones: 



KHa 



MI3 



3695360 
9 : 



96, 



l|2 
|u 4 || 3 



Ml 



6912, 
L366793248 



675 



1-^2 11 i 



I l|2 
P5|| 3 



45440, (5.10) 

2243123779689032 
186046875 ' 



11^*52 1| 3 is a ratio of integers where the numerator and the denominator have 19515 
and 19463 digits, respectively. Table 1 reports 1 1 it j 1 1 § for j = 0, ..., 52, in the 16 digits 
decimal representation. 

Table 1. The squared norms 1 1 1* ^ 1 1 3 - 



3 


ll?y-ll 2 
Il"jll3 


3 


ll?y-ll 2 
Il"yll3 


1 


6912 




27 


4.070323867244879 x 10 27 


2 


45440 




28 


3.800202819232687 x 10 28 


3 


4.105955555555556 x 


10 5 


29 


3.554589555246873 x 10 29 


4 


2.024878885925926 x 


10 6 


30 


3.330557264153261 x 10 30 


5 


1.205676676745595 x 


10 7 


31 


3.125627141295364 x 10 31 


6 


8.452219877103332 x 


10 7 


32 


2.937654907691943 x 10 32 


7 


6.152775603322622 x 


10 8 


33 


2.764771414352126 x 10 33 


8 


4.791192836997696 x 


10 9 


34 


2.605347861791808 x 10 34 


9 


3.628869598772102 x 


10 10 


35 


2.457968790658826 x 10 35 


10 


2.825486371143428 x 


10 11 


36 


2.321406184470901 x 10 36 


11 


2.228507964437443 x 


10 12 


37 


2.194593722846032 x 10 37 


12 


1.821213808657725 x 


10 13 


38 


2.076602420620089 x 10 38 


13 


1.539790191793044 x 


1Q 14 


39 


1.966618988613002 x 10 39 


14 


1.341372343677860 x 


10 15 


40 


1.863927582086700 x 10 40 


15 


1.190159209731028 x 


10 16 


41 


1.767894900465337 x 10 41 


16 


1.066432595016119 x 


10 17 


42 


1.677958174980847 x 10 42 


17 


9.598519025230687 x 


10 17 


43 


1.593615440091581 x 10 43 


18 


8.662788463495777 x 


10 18 


44 


1.514417532673484 x 10 44 


19 


7.840631870939454 x 


10 19 


45 


1.439961389630372 x 10 45 


20 


7.122921654632158 x 


10 20 


46 


1.369884345517744 x 10 46 


21 


6.499436510134908 x 


10 21 


47 


1.303859232337703 x 10 47 


22 


5.957837347113741 x 


10 22 


48 


1.241590147950303 x 10 48 


23 


5.485035371335649 x 


10 23 


49 


1.182808795435820 x 10 49 


24 


5.068929708200902 x 


10 24 


50 


1.127271314453561 x 10 50 


25 


4.699401376031744 x 


10 25 


51 


1.074755536205362 x 10 51 


26 


4.368534165204974 x 


10 26 


52 


1.025058601409640 x 10 52 
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Let us pass to the squared norms ||w^(£)||3- As an example, the result for iV = 5 is 

n (k\, v„9 9 258304 A 104566912 , 
« () (*) 3 = 96 + 6656£ 2 + — — £ 4 + — t 6 5.11 

9513575648 8 2243123779689032 in 

t 8 + 1 10 . 



70875 186046875 

There is no room to report here the results obtained for all the other values of N, 
especially in the rational form for the coefficients. However, we can write some of 
them in the 16 digits precision; in particular, 

\\u m (t)\\l (5.12) 

= 96 + 6656 1 2 + 2.870044444444444 x 10 4 t 4 + 1.993359937918871 x 10 5 1 6 

+1.058054454761424xl0 5 £ 8 +1.781444415306641xl0 6 £ 10 +2.740017914111055xl0 6 £ 12 

-7.321985472578865xl0 6 £ 14 +4.183410651491110xl0 6 £ 16 +1.457483700816015xl0 8 £ 18 

-1.768517246168822xl0 8 £ 20 +4.196205149715839xl0 8 £ 22 +3.648789154816725xl0 9 £ 24 

-2.178830191383206xl0 10 £ 26 -1.394064522752687xl0 10 £ 28 +2.954202883502504xl0 11 t 30 

+1.283692616423054x 10 11 t 32 -4.543575106022102x 10 12 1 34 +4.789569007452901 x 10 12 1 36 

+2.830635227431622xl0 13 t 38 +4.470168139346678xl0 13 t 40 -6.910532995061547xl0 14 t 42 

+1.457019276470951 x 10 14 t 44 +9.053007124662626x 10 15 1 46 -8.939780851014422 x 10 15 1 48 

-1.019952729404346xl0 17 t 50 +1.137772938577812xl0 17 t 52 +1.644161010427522xl0 18 t 54 

-4.571936581656874xl0 18 t 56 -3.140936865806385xl0 19 t 58 +2.408085513008218xl0 21 t 60 

-1.107900217253947xl0 23 t 62 +4.186064092726056xl0 24 t 64 -1.674853723772203xl0 26 t 66 

+6.911508987260593x 10 27 1 68 -2.698282390313396 x 10 29 1 70 +9.951375797771149 x 10 30 1 72 

-3.558771163372845xl0 32 t 74 +1.232107326257251xl0 34 t 76 -4.045044388392564xl0 35 t 78 

+1.242344004413423 x 10 37 1 80 -3.561397641466941 x 10 38 t 82 +9.520206481050174x 10 39 1 84 

-2.357932432543021 xl0 41 t 86 +5.354229494719748xl0 42 t 88 -1.103667607665446xl0 44 t 90 

+2.052382635232918xl0 45 t 92 -3.436006560519912xl0 46 t 94 +5.184487278969682xl0 47 t 96 

-7.072466985323957x 10 48 1 98 +8. 759614973466463 x 10 49 1 100 -9.896987665647683 x 10 50 1 102 

+1.025058601409640 x 10 52 t 104 . 

The rest of the paper reports a number of facts stemming from our computations, 
with the interpretation that we suggest for them. 
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Verification of the outcomes of [5 J on || ix.^ Jv ^ (t) ||g. Our computations based on 
the systematic use of rational numbers have given essentially the same results as in 



[5] about 



\u 



(AT) 



\U 



(TV) 



(0.32)|| 



as a function of N, in the two cases t = 0.32 and t = 0.35. So, 

| u (7V) (0.35) l^ grows 



3 seems to approach a limit value for large N, while 



rapidly with N; our use of rational coefficients ensures that such a rapid growth is 
not due to cumulative rounding errors. In Figures 1-2, we report \\u 



(TV) 



ClS db 



function of N G {0, ..., 52}, in the two cases t = 0.32 and t = 0.35; these figures are 
very similar to the ones at the bottom of pages 235 and 236 of [5], respectively (but 
comparison requires a rescaling, since the H 3 norm employed in [5] differs from ours 
by a constant factor). 

We agree with [5] in interpreting these results as indications that the power series 
for this initial datum has a finite if 3 -convergence radius r 3 , with r 3 G (0.32,0.35). 
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Figure 1. ||uW(0.32)||§ as a function of 
iVe{0,l,...,52}. 



Figure 2. ||uW(0.35)||§ as a function of 

N e {0,1,..., 52}. 
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Further evidence on the i/ 3 -convergence radius of the power series. This 
comes from the root test (I3.16P for n = 3: 



73 



liminf ||uj|| 3 

j—t+oo 



Vi 



(5.13) 



Figure 3 represents \\uj\\ 3 as a function of j G {1, ..., 52}. For j = 36, 38, ..., 52 we 
have a very good interpolation 



1 11-1/J 

\'i 1 -II ' r^j 

I u jll3 — 



0.32158 



1.20125\ 



1.38458 



(5.14) 



(obtained assuming for the interpolant the form A — (B/j) c , and applying the least 
squares criterion). The right-hand side of (15.14D approximates ||wj|| 3 



-Vi 



with a mean 



quadratic error < 10 5 (averaging, as indicated, for j = 36,38, ...,52; if we average 
over the larger range j = 16, 18, ..., 52, the mean quadratic error is < 10~ 4 ). 

Assuming that (j5.14p approximates HM7II3" with a similar precision for arbi- 



trarily large j, but keeping prudentially only two digits in our final estimate, we 
conclude with an estimate 

0.32 < r 3 < 0.33 . (5.15) 



0.35 




0.3 


.-- — "" 
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- /^ 
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- 
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0.05 


j 



10 
Figure 3. \\u 



j 113 



20 30 40 50 

1/3 as a function of N e {1, ..., 52}. 
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Reminder estimates for the series expansion of u(t) in H 3 . Let iV G 

{0, 1, 2, ...}; of course 



hOO 



U(t) -U iN \t) = J2 U ft fOT t G (" T 3 5 r 3) ; ( 5 - 16 ) 

j=N+l 



this implies 



+oo 



||«(t)-«W(t)|| 3 < J2 IMW fortG(-r 3 ,r 3 ) . (5.17) 

j=N+l 

To go on, we need a guess on the behavior of the norms 1 1 w^ 1 1 3 . To this purpose, let 
us consider the sequence 

fjL 3j := 0.32^' || Uj || 3 (j = 0,1,2...), (5.18) 

recalling that 0.32 is the lower bound for r 3 in (j5.15p . From the norms available up 
to j = 52, we can check that (/%) is decreasing while j ranges in {1,3, ...,52}; by 
extrapolation, let us assume that (/%) is decreasing on the infinite set {1,2, ....}. 
So, /i 3 j ^ /x 3 at for integer j ' ^ N ^ 1, i.e., 

Ik,- lis ^ -^ for j ^ iV ^ 1 . (5.19) 

" J " 0.32^ 

For £ G (—0.32, 0.32), inserting this inequality into ( I5.17P we get ||w(t) — w^(t)||3 ^ 
/^siv E& + i IV0.32p' = ^ 3N \t/0.32\ N ^ ES IV0.32p', i.e., 

It/0 32l 7V+1 
H^)-^W(t)|| 3 ^^3iv }_; ' 32| for tG (-0.32,0.32), iVG {1,2,3,...}. (5.20) 

Of course, this is a conjecture based on the previous extrapolation. For the practical 
application of the reminder estimate (15.201) . we mention that (rounding up from 
above) 

/i 35 = 11.7, u, 3W = 5.99, /i 320 = 3.39, (5.21) 

/i 330 = 2.61, /! 34 o = 2.20, yu 352 = 1.88 . 

No blow-up at Ts- After accumulating indications that the Taylor series for u(t) 
has an if 3 - convergence radius r 3 G (0.32,0.33), in the rest of the section we will 
present evidence that u(t) does not blow up at t = r 3 . 

The power series for ||w(t)||g; an indication that u(t) exists up to t = 0.47 
at least. The results (13.191) ( 13.201) with n = 3 give a formal series expansion 

+oo +oo j 

II 5Z u ^ 3 Us = 5Z " 3 ^' U3 i := y^(^kj-^3 e R, ^3j = for j odd ; (5.22) 
i=o j=o e=o 
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the series X^^o v ?>ft ^ as a convergence radius 



^ 3 = liminf |z/ 3i |- 1/J . (5.23) 

Recalling that (—T,T) is the domain of the solution u, we know (from (I3.23P ) that 

+00 
t 3 ^9 3 <:T, \\u(t)\\i = Y,v^ for te(-6 3 ,9 3 ). (5.24) 

3=0 

In the sequel, for N = 0, 1, 2, .... we also consider the partial sums 

N 

(*):=£>,■*»'. (5.25) 

j=0 



,(N) 



Of course, «W(t) = J2f=o u i tJ is such that «(*) = « W ( f ) + 0(t iV+1 ) for t -> 0; this 
implies ||u(t)||! = \\u^{t)\\l + 0(t N+l ), whence 

i/$ JV) (t) = ||u (JV) (t)||l , c , Ar • (5.26) 

3 w " w " 3 t fe ^0 for k > N v ; 

With this remark, the previuos computations of ||u^(£)||3 up to iV = 52 also give 
the partial sums z/ 3 (£) for N = 1, ...,52 or, equivalently, the coefficients u 3 j for 
j = 0, ..., 52. For example, 

ofi rrer 258304 2825587712 

z/ 30 = 96, ^32 = 6656, z/34 = — - — , z/ 36 = — — — — , (5.27) 

9 14175 

52545219363488 10025320340466597351685768 

^38 — TT^TT-, lr >r ' ^310 



496621125 5627635784943046875 

1/352 is a ratio of integers where the numerator and the denominator have 2610 and 
2593 digits, respectively. 

The 16-digits representation of the coefficients v 3 j for all j £ {0, ...,52} can be 
obtained from Eqs. ( 15.1 2p ( 15. 26ft ; more precisely, 

u 3j = coefficient of P in (15TT21 . for j = 0, ..., 52 . (5.28) 
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From the above data, one can try to make predictions on the convergence radius 
#3 of the series X^j^o u 3j^ ■ m Figures 4-7 we report the partial sums v\ (£) as 
functions of N e {0, ..., 52}, in the four cases t = 0.45, 0.50, 0.55, 0.60. For t = 0.45, 
the function N \— > u^ ' (t) seems to approach a limit value for large N. The situation 
is not clear for t = 0.50, due to the appearing of small oscillations; for t = 0.55 and 
t = 0.60, the oscillations of N \— >■ v\ (t) are large and their amplitude increases with 
N. We regard these results as indicating that XlS> v ?>ft ^ s convergent for t ^ 0.45 
and not convergent for t ^ 0.55; in other words, for the convergence radius we have 
a conjectural estimate 

0.45 < 3 < 0.55 . (5.29) 




30 
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Figure 4 



(0.45) as a function of 




N e {0,2, ...,50, 52}. 



Figure 5. ^3 (0.50) as a function of 
N e {0,2,..., 50, 52}. 




Figure 6 




(0.55) as a function of 



TV e {0,2, ...,50, 52}. 
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Figure 7. 1/3 (0.60) as a function of 
N e {0,2,..., 50, 52}. 
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Another way to estimate #3 comes from the root test (j5.23p . Figure 8 is a graph 
of |^3j| -1 ^ as a function of j G {2,4, ..., 50, 52}. For j = 36, 38, ..., 52, there is a 
fairly good interpolation 



m\ 



- 1/j ~ 0.484 



8.4J 



2.19 



(5.30) 



(obtained assuming for the interpolant the form A — (B/j) c , and applying the least 
squares criterion); here, the right-hand side approximates |^3j| -1 ^ with a mean 
quadratic error < 0.01 (let us repeat it, for j between 36 and 52). Assuming that 
the above interpolant behaves similarly for all larger (even) j, and considering 6* 3 = 
lim inf j_>. +00 |^ 3 j| _1 / J we are led to use 0.484 ± 0.01 as upper and lower bounds for 
it; rounding up to two digits we obtain the inequality 

0.47 <9 3 < 0.50 , (5.31) 

which is compatible with f l5.29[) . Now, recalling that #3 is a lower bound on the time 
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Figure 8. l z/ jl3 ±/J a s a function of j <G 
{2,4, ..., 50, 52}. The dashed line is the graph of 
the interpolant in (|5.30p . for j E [36, 52]. 



of existence T of the solution u (see (j5.24p ). we are led to the final estimate 

0.47 < T SC +00 . (5.32) 

In particular, as anticipated, we have indications that u does not blow up near the 
if 3 -convergence radius r 3 . 
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6 Possible blow-up at larger times for the Behr- 
Necas-Wu datum, via Pade approximants 

A few words on Pade approximants. Let us be given an analytic function 
/ : / — y C, t !->■ f{t), with / a neighborhood of zero in Ror C. Let p,q <E {0, 1, 2, ...}; 
we recall that the Pade approximant of order (p, q) of /, if it exists, is the unique 
complex function \p/q]f = [p/q] of the form 

ip/q](t) = a o + M + - + V P , (6.1) 

such that 

f(t) = \p/q](t) + 0{t p+q+l ) for t -> ; (6.2) 

the above condition determines the p + q + 1 unknown coefficients oo, ..., b q as func- 
tions of the derivatives /^(0), j = 0, ...,p + g; the domain of [p/q] is the largest 
subset of C where the above ratio is defined. The family of all approximants [p/q] 
(j),q = 0, 1, 2, ...) forms the so-called Pade table of /; the approximants with p = q 
are called diagonal. 

There are several results and conjectures about the convergence to / of the 
Pade approximants [p/q] with p or p, q large. In particular, the so-called "Pade con- 
jecture" (or "Baker- Gammel- Wills conjecture") states that, for a meromorphic func- 
tion / on a disk of C, there is a subsequence [pi/pi] (£ = 1,2,3,...) of diagonal 
Pade approximants that, for £ — y +oo, converges to / uniformly on each compact 
subset of the disk minus the poles of /. This conjecture has been proved for special 
classes of meromorphic functions (see jl] [21] [22] and references therein). 

It is found experimentally that the Pade approximants of large order (and, in 
particular, the diagonal approximants [p/p]) work as well for many non meromorphic 
functions, describing accurately their behavior even close to non polar singularities. 

Pade approximants for ||w(t)||g, and possible evidence for a blow-up. The 

previous considerations can be applied (for suitable n) to the function /(£) := 
||w(£)||^, where u is the solution of the Euler equation with a given datum uq. 

One can ascribe to a number of works the idea of using the Pade approximants 
for such a function; as in the Introduction, we mention [5] [12] [U] [2D] (and some 
references therein). As already remarked, these papers have considered initial data 
uq different from the one of Behr-Necas-Wu (e.g., the Taylor-Green vortex); further- 
more, they have generally considered the Sobolev norm of order n = 1. 

Here we are focusing on the (maximal A-) solution u for the Behr-Necas-Wu 
datum; from now on, [p/q] stands for the Pade approximants of the analytic function 

t^f(t):=\\u(t)\\j. (6.3) 
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We conjecture that, for certain large p, [p/p] approximates well the function t i— >■ 
\\u(t) || 3 (and even its analytic continuation to the complex plane). From the previous 
paragraphs, we have the derivatives /^(0) = j\ v^j for J ' — 0, ..., 52; this information 
suffices to determine all the Pade approximants [p/q] for p+q ^ 52 and, in particular, 
all the diagonal approximants [p/p] for p — 0, 1, ..., 26. 

It turns out that the diagonal approximants [p/p] exist in the cases of even 
order p = 0,2, ...,26, while they do not exist in the odd cases p — 1,3, ...,25 (the 
reason being, essentially, that the power series for f(t) about zero contains only even 
powers oft). Let us consider, for example, the approximant [12/12]. Its numerator 
and denominator are polynomials with rational coefficients, too large to be written 
explicitly; however, we can use the 16-digits approximation for the coefficients and 
write 

[12/12IW = gg , (6.4) 

N 12 (t) := 96 + 6.680481407149543 x 10 3 1 2 + 3.08095009988031 x 10 4 t 4 

+2.3462351635051233 x 10 5 £ 6 + 2.407391215430808 x 10 5 1 8 

+2.5575522886490226 x 10 6 t 10 + 3.094974424148063 x 10 6 t 12 , 

D 12 (t) := 1 + 0.255014657807743 t 2 + 4.288322833232482 t A - 5.985294148961588 1 6 

+8.973150435320479 1 8 + 66.29326162173366 1 10 - 612.1107629833056 t 12 . 

The poles of [12/12], which are the zeros of Di 2 , are simple and occur at the points 
t = ±0.294020 ± 0.464361 i (\t\ = 0.549617) ; (6.5) 

t = ±0.511609 ± 0.301416 i (\t\ = 0.593797) ; 
t = ±0.606004 i , t = ±0.626199 

(here and in the sequel, ± means that we can choose independently the signs for 
the real and imaginary part, e.g., + for the real and — for the imaginary part). So, 
the singularities of minimum modulus of the approximant [12/12] are at anyone of 
the points T = ±0.294020 ± 0.464361 i, such that \T Q \ = 0.549617; furthermore, the 
real singularities closest to the origin are at anyone of the points T* = ±0.626199. 

We have performed a similar analysis for all the approximants [p/p], with p = 
14, 16 26; the results are summarized in Table 2. 
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Table 2. Poles of the Pade approximants \p/p](t) to ||ii(£)||j 



T 
T, 



pole closest to the origin (with modulus |T |); 
real (or almost real) pole closest to the origin. 



\p/p] 


T 

- 1 


\T 1 
1 ° 1 


T 


[12/12] 


±0.294020 ± 0.464361 i 


0.549617 


±0.626199 


[14/14] 


±0.281333 ± 0.445002 i 


0.526474 


±0.656185 


[16/16] 


±0.283300 ± 0.446498 i 


0.528790 


±0.661087 


[18/18] 


±0.283081 ± 0.445859z, 


0.528134 


±0.660118 


[20/20] 


±0.345307 ± 0.348713 i 


0.490752 


±0.621387 ± 0.047708 i 


[22/22] 


±0.350239 ± 0.350695 i 


0.495635 


±0.541967 


[24/24] 


±0.349063 ± 0.350777 i 


0.494863 


±0.609804 ± 0.0383530 i 


[26/26] 


±0.0714399 ± 0.508700 i 


0.513692 


±0.816133 



Let us point out some features of the Pade approximants Table 2, with their possible 
implications: 

(i) For all the approximants [p/p] in the table, the poles of minimum modulus 
occur at points T with |T | ~ 0.5. There is not a clear trend of |T | as a function 
of p, so we limit ourself to consider the mean of |T | for p = 12, ...,26 which is 
(|T |) = 0.515995, with a mean quadratic error AT < 0.02. On the other hand, for 
a holomorphic function, the convergence radius of the power series centered at zero 
is the modulus of the singularity closest to the origin. So, assuming that the above 
[p/p] describe approximately the singularities of fit) = \\u(t)\\^, we can derive from 
these approximants an estimate of a convergence radius #3 for the power series of 
f(t). More precisely, assuming |T | — AT < 9 3 < \T \ ± AT and rounding up to 
two digits, we obtain from the above Pade approximants an estimate 



0.49 < 0, < 0.54 



(6.6) 



this is compatible with the estimate on 83 obtained in Section by other means (see 
Eq. (15.311) and the discussion before it). 

(ii) The [p/p] approximants of Table 2 have real poles (symmetric with respect 
to the origin), with the exceptions of [20/20] and [24,24] which, however, possess 
"almost real" poles, close to the real axis (P)- In the table, we have denoted with 
T* the real (or almost real) singularities closest to the origin. The mean of [T*! for 
for p = 12, ...,26 is (|T*|) = 0.649489, with a mean quadratic error AT* < 0.08 
(however, there are large deviations from the mean in the special cases p = 22 and 
p = 2Q). 



7 The occurring of almost real singularities has also been pointed out in [T2] [30] while analyzing 
the Pade approximants for ||u(i)||i, with initial conditions uq different from the Behr-Necas-Wu 
datum. 
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The above results on the singularities T* somehow suggest that /(£) = ||w(£)||§ could 
diverge for t — > T~ (and t — > (— T) + ), for a suitable T; if we assume for T the upper 
and lower bounds |T*| ± AT*, rounding up to two digits we get 

0.56 < T < 0.73 . (6.7) 

If such a conjectured divergence of f(t) actually occurred, the solution u of the Euler 
equation with the Behr-Necas-Wu datum would blow up at T (and — T); admittedly, 
the indications for such a blow up are very weak. 

D-log Pade approximants for ||w(t)||g. As well known, the D-log Pade ap- 
proximants of a function t (->■ f(t) are the Pade approximants for the logarithmic 
derivative /// (' := d/dt). These approximants are generally regarded as more suit- 
able for describing the behavior of / close to singularities, even of non polar type. 
In particular, the presence of a singularity at a point T*, say real, and a behavior 
of the type [p/p]r,r ~ A*/(T* — t) for t — > T~ is regarded as an indication that 
f(t) ~ const/ (T - t) x for real t^T~, where T ~ T* and A ~ A, pQ. 

We have attempted an analysis of the function f(t) := \\u(t)\\l via the approx- 
imants \p/p]f/f, with odd p ^ 25 |j); the results are very unstable with respect to 
the order, and ultimately not sufficient to get any indication of blow-up. (p) 



8 Our function has the form f(t) = F(t 2 ); in such a case, for odd p, the D-log approximant of 
/ of order (p,p) is (up to a factor 2t) the D-log approximant of the function s H> F(s) of order 
(p/2 — 1/2, p/2 — 1/2). On the contrary, for even p, the (p,p) D-log approximant of / cannot be 
interpreted in terms of F. 

9 Here is a more precise description of the computational outcomes. The D-log approximants 
of order (p,p) for p = 17,19,21 have real singularities at points T» ~ 0.72 and are such that 
[p/p\f/f ~ ^*/{T* — t) for t h-> T~, with A» ~ 2.6; so, for ||u(t)|| 3 = y/f(t) we have a conjecture 
||u(t)|| 3 - const. /(T - t) a with T ~ 0.72 and a ~ A*/2 ~ 1.3. This value of a agrees with the 
Bcale-Kato-Majda bound a ^ 1 in the event of blow-up (see Eq. (|2.16p ); it agrees as well with the 
(conjectural) bound a ^ 6/5, obtained extrapolating from R 3 to T 3 the estimate (12.171) . 

On the contrary, the D-log approximant of order (23, 23) for / has no real (nor almost real) 
singularity. Finally, at the order (25, 25) there is a real singularity for T, ~ 0.52, and [25/25] t ,* ~ 
A*/(T„ — t) for t M> T~, with —0.002 < A* < 0.002 (there are numerical difficulties in a more 
precise determination of A»). Returning to ||u(i)||3 = \/f(t), the [25,25] Pade would suggest 
\\u(t)\\ 3 - const. /(T - t) a with T ~ 0.52 and -0.001 < a < 0.001. This statement is an absurdity 
even in the case < a < 0.001, since it contradicts the Beale-Kato-Majda bound (J2.16I) a ^ 1. 
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7 Conclusions 

The previous results about the Behr-Necas-Wu datum uq support our statements in 
the Introduction, i.e.: 

(a) The power series for u has an H 3 convergence radius r 3 such that 0.32 < r 3 < 
0.33 (see Eq. ( jg35j) ). 

(b) There is no blow-up at time r 3 and the (maximal A-) solution u of the Euler 
Cauchy problem exists, at least, up to a time (9 3 (the convergence radius for the 
series expansion of ||u(£)|| 3 ), for which we have from f !5.3ip the estimate 8 3 > 0.47. 

(c) The Pade approximants for ||w(£)|| 3 in Table 2 give weak indications that u might 
blow up at a time T, with 0.56 < T < 0.73 (see Eq. (16. 7p ). 

We think that the evidence given in this paper is rather strong for (a)(b). As 
for (c), doubts on the blow-up conjecture arise not only from the rather erratic 
behavior of the real singularities in the computed Pade approximants; in fact there 
are more general reasons, recalled at the end of the Introduction, suggesting caution 
in deriving blow-up results from the Pade approximants. 
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